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Abstract 

In a solvable model of two dimensional SU (N) (N — > oo) gauge fields interacting with 
matter in both adjoint and fundamental representations we investigate the nature of the 
phase transition separating the strong and weak coupling regions of the phase diagram. 
By interpreting the large N solution of the model in terms of SU(N) representations it is 
shown that the strong coupling phase corresponds to a region where a gap occurs in the 
spectrum of irreducible representations. We identify a gauge invariant order parameter 
for the generalized confinement-deconfinement transition and give a physical meaning 
to each phase in terms of the interaction of a pair of test charges. 
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1 Introduction 



The study of systems where interactions are mediated by non-Abelian gauge fields is of direct 
relevance to the physically interesting case of quantum chromodynamics (QCD). At high 
temperature or density these systems are expected to undergo a phase transition where the 
character of the effective degrees of freedom changes dramatically. For example, in the low 
temperature phase of four dimensional QCD, quarks and gluons carrying colour charge are 
not observed but rather confined into composite baryons and mesons. It is expected, and can 
be shown in numerical simulations on the lattice, that at sufficiently high temperatures this 
confinement is relaxed and the fundamental degrees of freedom become mobile in a quark- 
gluon plasma. Quantifying the differences between these phases has been a subject of study 
for some time now an d is adequately understood only in the case of pure gluo-dynamics 
without quarks. Here the Polyakov loop operator [Q] 

< Tr g(x) >=< Tr P rf^ M^)dr > ^ ^ 

provides an effective order parameter H for the transition from the confined to the deconfined 
phase by testing to see if the symmetries of the action are realized faithfully. As we will 
show it is useful to consider the trace of the group element g in group theoretic terms as 
defining a group character. Taking g in different irreducible representations will allow us to 
unambiguously define the strong and weak coupling regimes of a two dimensional model even 
in the presence of fundamental matter. 

Solvable models are often of use in developing new ideas and testing hypotheses. Here 
we will use the model of the non-Abelian Coulomb gas in one space and one time dimen- 
sion which has a non-trivial phase structure qualitatively similar to what we expect in four 
dimensional QCD. In particular we will discuss a finite temperature model of static adjoint 
and fundamental representation charges interacting through SU(N) gauge fields in the limit 
N — > oo 0. It is most conveniently presented as a unitary matrix model and its solution 
can be found explicitly in the large N limit by calculating the element g of the gauge group 
which dominates a saddle-point approximation of the partition function. The phase structure 
of the model, as we will show, consists of two general regions, one which we will define as 
weak coupling where go ~ 1 and one of strong coupling with go away from the identity ele- 
ment. It was noted [|J] that the behaviour in the system of Tr g$ is markedly different in each 
phase when varying the parameter n. Some effort was made to interpret these characteristics 
physically and using the language of group characters we complete that task here. 

Even though two dimensional Yang-Mills is a dynamically trivial theory, as the rank of the 
symmetry group is taken to infinity, group theory can drive phase transitions. Transitions of 
this type were first noted long ago in the lattice theory |J but these are now considered to be 
lattice artifacts. More recently such phase transitions have been noted in the continuum with 
the Douglas-Kazakov transition || on the sphere and the related transition on the cylinder 
J7[ being prime examples. In these cases the theory is solved for large rank symmetry group 
in terms of a single irreducible representation which saturates the evaluation of the partition 
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function in a saddle-point approximation. The phase transition corresponds to a point where 
the distribution of occupation numbers for the rows of the associated Young table develops a 
gap B0. 

In the present case under consideration the situation is somewhat different. The saddle- 
point is not determined in general by a single irreducible representation of the gauge group but 
by a linear combination of irreducible representations. This feature is also shared by Abelian 
and non-Abelian Coulomb gases in two dimensions with U(N) and SU(N) finite rank gauge 
groups H and can be generalized to the case of any compact Lie gauge group. In each of these 
cases the state vector of the system, \l/ is a class function and therefore can be represented 
by a linear combination of characters, xr °f the irreducible representations, R of the gauge 
group with coefficients an that depend on the parameters of the model (temperature, pressure, 
gauge coupling constant...) 

ng(x)}=J2dRXR(9(x)) (1.2) 

R 

Consequently we see that there are two different points of view to take with solving these 
models in two dimensions. One is to find a dominant configuration of the gauge group go 
and the other is to find the dominant linear combination of irreducible representations, \&. 
The main objective here is to quantify the connection between these two views and use it 
to characterize the differences between the strong and weak coupling regimes of the non- 
Abelian Coulomb gas. As we will see, the characters of the gauge group are completely 
determined by traces of powers of the gauge matrices, Tr g n . In the two dimensional model 
under consideration we will show that the vanishing of particular coefficients an provides a 
convenient way to characterize the different phases of the model. Clearly, if a particular is 
vanishing then the system does not have excitations which can effectively screen a charge in 
irreducible representation R interacting with its conjugate R. In this way we will be able to 
identify an order parameter for the transition from strong to weak coupling and give a physical 
definition of the confinement-deconfinement transition with fundamental matter present. 

The layout of this paper is as follows: First we present a short description of the non- 
Abelian Coulomb gas model which will be used as a test-bed for our program of using group 
characters to describe the different phases of a system of interacting gauge fields. In particular 
we will show that this model possesses an interesting phase structure. We follow with an 
explicit demonstration of the connection between a gauge group element g and the characters 
of the irreducible representations of the the gauge group. Applying these general results to 
the case of the model at hand we will show that the spectrum of irreducible representations 
provides a clear way to distinguish between phases of the model with fundamental matter in 
much the same way the Polyakov loop operator does in pure gluo-dynamics. 

2 Review of the Non-Abelian Coulomb Gas 

The model on which we will base our investigation is that of the non-Abelian Coulomb gas 
in two dimensions as previously studied in This model, while lacking both dynamical 
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gauge and matter degrees of freedom, exhibits many of the features expected || in higher 
dimensional systems. For instance in SU(N) gluo-dynamics for N > 3 we expect a first 
order deconfmement transition with or without adjoint matter. As one adds matter in the 
fundamental representation it is expected that the latent heat associated with this transition 
drops until the transition becomes continuous or is completely washed out. As we will show 
in this section, the two dimensional non-Abelian Coulomb gas shares this behaviour and is 
explicitly solvable in the large N limit. The fact that we consider the limit where the rank 
of the symmetry group is taken to infinity is necessary to generate a phase transition in this 
simple model. 

We will now outline the details of this model which are relevant to the phase structure 
and our interpretation of it. Beginning with the canonical quantization of 1 + 1 dimensional 



Yang-Mills theory with gauge coupling e at finite temperature T it can be shown |L0|, [11], [| 
that the grand partition function for a system of interacting static colour charges is that of a 
gauged principal chiral model 

Z[X,T] = J [dA][dg] exp - J dx Tr \Vg + i[A,g]\ 2 - A |Tr g\ 2 - 2Nk ReTr g 

7 y (2.1) 
Here the integration over gauge fields A effectively enforces Gauss' law as one integrates over 
all elements of the gauge group with the Haar measure [dg]. The fugacities of the adjoint and 
fundamental charges are given by the parameters A and Nn, respectively. Since we consider 
the matrix- valued fields A and g to be taken in the fundamental representation of SU(N), the 
large N limit will lead directly to the familiar situation of matrix models with large N x N 
matrices. In order to keep all terms in the action of ( |2.1| ) at leading, N 2 , order in this limit 
we will restrict parameters of the system such that 7 = J^, A and k are each of 0(1). 

For the discussion of a confinement- deconfmement phase transition the most important 
aspect of the action in ( |2.1| ) is a global symmetry S[A, g] = S[A, z g] when the fundamental 
charge fugacity k vanishes. Here z is a constant element from the center of the gauge group, 
which for U(N) is U(l) and for SU(N) is Zjy. It is this symmetry and its (thermo-) dynamical 
breaking that leads to the deconfmement phase transition in this model. If k 7^ the question 
of what remnants of this symmetry persist is one we will answer in the next sections. 

Additionally, there is a gauge invariance that can be used to diagonalize the matrices 
gij(x) = e lai ( x '5ij. The density of eigenvalues p(9,x) = jfJ2iLi$(8 — cti(x)) corresponding 
to the large TV saddle-point evaluation of ( |2.1| ) now completely characterizes the properties of 
the system. Our goal is to find this distribution of eigenvalues. Without loss of generality we 
can consider the Fourier expansion 

p(B,x) = ±- + ±- £ c n (x)e- m9 , c n (xy = c. n (x) (2.2) 

ZTT Z7T ~f„ 



The configurations of the eigenvalue density (|2.2j) that saturate ( |2.1|) at large N can be found 
via the collective field theory approach [13|, |lOj. The method is essentially based on the 



relation between matrix quantum mechanics and non-relativistic fermions [12]. Leaving the 
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details to [|J, it can be shown that a solution of the saddle-point evaluation of ( |2.1| ) is given 
by 



(0) = I V V E + 2 ^ Cl + K ^ cos0 where pis real 
PoU \ otherwise 1 ' ' 

The constant of integration E has physical interpretation as the Fermi energy of a collection 



of N fermions [|12j on the circle subject to a periodic potential. It is fixed by requiring the 
eigenvalue distribution to have unit normalization. Furthermore, since the potential due to 
adjoint charges is non-local in eigenvalue space, the Fourier moment C\ (see (|2.2|)) must be 



self-consistently determined [14 



ci = dO p o (0)cos0 (2.4) 



This pair of conditions is most conveniently analyzed by introducing a new parameter fi = 
E/{2(\c\ + k)) and the integrals over the positive support of /i + cos 9, 

I n (fi) = — J dO cosnO^fi + cosO (2.5) 
In terms of /i, the solution of the normalization and moment conditions is given by 

d = W\ (2.6) 



and 

k 1 A Ii (/J,) 



(2.7) 



7 AIM 2 

This last relation gives a family of lines in the (A/7, k/t) plane parameterized by p. As 
is shown in Fig.l this family overlaps itself for lower densities of fundamental charges, k/j 
signalling the fact that there are multiple solutions to the equations of motion in this region 
of the phase diagram. Considering the free energy one can show Q that for lower densities 
of adjoint charges the stable solution has fi > 1 while at higher adjoint densities the stable 
solution has fi < 1. In the intermediate regime there is a first order phase transition. As the 
density of fundamental charges is increased the first order transition is smoothed out and a 
third order phase transition persists along the line — 

The parameter fi is now seen to be useful for two different reasons. First it characterizes 
the general structure of the phase diagram (Fig. 2) where the 'strong coupling' regime is the 
region with fi > 1 and the 'weak coupling' regime has /x < 1. As well, and of more importance 
for our analysis, we find that the expectation values of traces of powers of the group element 
g are given as a function of the single parameter \i 

< Tr g»/N >= c n = ^ (2.8) 

Consequently, it makes sense for our purposes to re-define the eigenvalue distribution in terms 

Of fjL 



po(0,aO = —rr^v^ + cose ( 2 - 9 ) 

In the next section we will use this definition and its connection to the dominant configuration 
of the gauge element to analyze the phase diagram in terms of group theory. 
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Figure 1: Plot of the lines (2.7) for /j, ranging from 0.4 (upper right corner) to 75 (line at 
the extreme left). The region of overlapping lines corresponds to a region of first order phase 
transition. 




K/y 



Figure 2: Schematic picture of the phase diagram. The doted curve marks the first order part 
of the critical line. The solid curves above and below it are the boundaries of the area with 
two possible phases. They join at a point which shows second order behaviour. For larger fc/7 
we find a third order line (fi = l) marked by a solid line. 
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3 Order parameters and some group theory 



As is known, in the case of pure gluo-dynamics, the realization of the center symmetry of the 
gauge group governs confinement @. The Polyakov loop operator Tr g(x), which is related to 
the free energy -Tlog< Trg(x)Trgt(0) > of a conjugate pair of static, external fundamental 
charges separated by a distance x, serves as an order parameter || to test confinement. Since 
Tr g(x) transforms under the center as Tr g{x) — > z Tr g(x), the expectation value of the 
Polyakov loop operator must average to zero if the center symmetry is preserved. Physically 
this suggests that an infinite amount of energy is required to introduce a single fundamental 
test charge into the system. The presence of a gas of fundamental charges (k ^ 0) changes 
this situation though by explicitly breaking the center symmetry. Consequently we lose the 
Polyakov loop operator as an order parameter for phase transitions in the system. In this 
section we introduce a suitable generalization of the Polyakov loop operator which will allow 
us to identify a new order parameter. 

As seen in the previous section, the solution of the non-Abelian Coulomb gas with adjoint 
and fundamental representation charges is completely characterized by a Fourier sum of the 
traces c n =< Tr g n /N > - the higher winding Polyakov loops. As noted in [Q the character 
of these traces changes between the strong and weak coupling regimes. In particular, in the 
strong coupling (// > 1) phase c n is damped exponentially with n while in the weak coupling 
(/i < 1) phase the damping follows a power law behaviour. Now we will reconsider this 
behaviour in terms of group theory. 

Since the matrix g is an element of the special unitary group, its trace in an irreducible 
representation, R defines the group character for that representation 

XR(g)= Tr R9 (3-1) 

For the N dimensional fundamental representation of SU(N), F, the group character is just 
the Polyakov loop operator described above since we are considering group elements to be 
taken in the lowest fundamental representation 

XF(g) = Trg (3.2) 

Further simple examples are the symmetric (S) and anti-symmetric (A) combinations of a 
pair of fundamentals where we have 

X s(g) = ^[(Trg) 2 + Trg 2 } , XA {g) = ~[(Tr gf - Tr g 2 } (3.3) 

A general relation between characters and the group elements is given by the Weyl formula but 
is not necessary for the following. A complete discussion can be found in standard references 
(see [[18|] for example). 

The main idea is that the eigenvalues of the group matrices, which are the only relevant 
dynamical variables in the grand partition function ( |2.1| ), are completely determined by the 
N quantities {Tr g n }, n = 1 . . . N. In turn these traces form an algebraic basis equivalent to 
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the characters of the N fundamental (completely anti-symmetric) irreducible representations 
of SU(N) (including the trivial representation). Here we will explicitly demonstrate the 
relationship between the basis of traces and the basis of group characters. Ultimately it is 
the group theoretic variables which we will use to characterize the phases of the model ( f2.1| ). 

The standard basis for general functions (of finite degree) of the eigenvalues of a matrix 
is the set of elementary symmetric functions {a r }. In terms of the eigenvalues Xj = e*^ of 
the group element g they are given by 



a 3 



j 

j<k 

j<k<l 



(3.4) 



a N 



Y[ Xj = det g — 1 



(3.5) 



with a r = for r > N. The relationship of the symmetric functions {a r } to the traces of the 
group elements, S n = Tr g n , is given |L5| by the determinant 



1 

jfc! 



Si 

s 3 

Sk-i 



1 

Si 
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Sk-2 

Sk-i 





2 

Si 

Sk-3 
Sk-2 





3 



s 2 

S3 





51 k-l 

52 S\ 



(3.6) 



Most importantly, it can be shown that the elementary symmetric functions are nothing 



more than the characters of the fundamental representations for the unitary group [15, 16 



That is, for the fundamental representation which is the anti-symmetric combination of k, N 
dimensional representations, Xk(g) = a k- 

The determinant fl3.6p can be evaluated ]T7[ in terms of a multinomial expansion most 
compactly stated in terms of a generating function 

<-l) k d k 



Xk(g) 



k\ dz k 



exp [- £ ™fz n ) 



n=l 



n 



(3.7) 



2=0 



For our purposes though it is useful to convert to a contour integral about the origin. 

-l) k r dz 



Xk(g) 



2ni 



yk + l 



exp 



71=1 



n 



(3.8) 



These last two expressions explicitly demonstrate the relationship between the group element 
g and the k th fundamental representation of the gauge group and are completely general 
results. 
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With these relations we see that there is a direct connection between the gauge group 
element g and the irreducible (fundamental) representations of the gauge group. In particular, 
in the previous section we have seen that in the large N solution of the non-Abelian Coulomb 
gas a certain configuration of the gauge matrix, g saturated the evaluation of the partition 
function (|2.1| ). Now it is natural to ask what is the configuration of irreducible representations 
corresponding to the dominant g . This corresponds to evaluating the expectation < Xk(g) > 
in the background of the non-Abelian gas. In principle this involves calculating expectations 
of the form < Tr g ni • • • Tr g Hr > but because of the factorization of gauge invariant objects 
in the limit iV — > oo, this reduces to a product of expectations, < Tr g ni > • ■ ■ < Tr g nr >. 
Consequently < Xk(g) > is determined by replacing Tr g n by its expectation value in (|3.8|) . Of 
course expectation values of the group element traces are intimately related to the eigenvalue 
density p(9,/i) (see (|2.2j)) hence, after performing an infinite sum, we obtain 



< Xa > \p{0,v)] = 



2ni 



f dz 


~N r 


j> — exp 


y J 







1 + z 2 - 2zcos9" 



y 2a 



(3.9) 



Note that we have defined a new real parameter a = k/N on the unit interval that effectively 
labels the fundamental representations in the large iV limit. Of course (|3.9|) now depends 
on a continuous variable and is of a slightly different functional form than the discrete case 
< Xk >■ In the remainder of this discussion we will consider only the character parameterized 
by a as defined in ( |3.9| ). 



4 Calculation of the expectation of fundamentals 

In this section we will concentrate on calculating < Xa > with eigenvalue density (|2.9|) for the 
non-Abelian Coulomb gas. This calculation will give a clear picture of the group theoretic 
excitations present in different regions of the phase diagram and consequently allow us to 
define an order parameter for the deconfinement transition even in the presence of fundamental 
matter. 

Since explicit evaluation of ( |3.9|) is difficult we begin with some special limiting cases. As 
\i — > — 1 the support of the eigenvalue distribution (|2.9|) vanishes at 9 = 0. The distribu- 
tion does not vanish though as it retains unit normalization and effectively becomes a delta 
function, 5(9). Consequently we find the gauge matrix g is just the identity at /i = —1, hence 

< Xa >= lim ( *)= 2 N \[J^ e- 2 ^- 1 ^ 2 (4.1) 

In this limit we find that the distribution of characters is symmetric about a = 1/2 as one 
would expect in a system where the total colour charge is vanishing. As well in this limit 
< Xa > is non-vanishing and all fundamental representations are present in the large N 
background solution of the model. As we will see, this result is generic in the weak coupling 
phase /i < 1. 
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In the opposite limit, as p — > oo, it can be shown that the eigenvalue distribution ( |2.9| ) 
approaches a constant value p = l/2n with the eigenvalues of the group element g becoming 
uniformly distributed on the unit circle. Since expectation values of the traces of powers of 
the gauge matrix are essentially Fourier transforms of the eigenvalue distribution, it is easy 
to see that < Tr g n >^ in this limit and 

< Xa >^ <V ( 4 - 2 ) 

This limit corresponds to the extreme strong coupling phase of the model where the Polyakov 
loop operator (< Tr g >~< Xi/n >) has vanishing expectation value and the standard 
analysis would point to a phase where colour charges are strictly confined into hadron-like 
structures. 



In general the integral (|3.9|) can be evaluated by saddle-point methods in the large N limit 
in which we are interested. The relevant action in this limit is 

S(a, fi,z)=J d6p(6, p) log ( — j (4.3) 

Solving the stationarity condition, dS/dz\ Z0 = 0, for a in terms of z we find the saddle-point 
condition for the large N behaviour of the integral (|3.9|) is given by the relationship 



f aq (a \ z o(z -cos9) 

a = / d6p{6,p)—— — 4.4 

J 1 + Zq — 2z cos 6 

Since a is a real parameter restricted to the unit interval [0, 1] it can be shown that the 
saddle-point value of the parameter z$ is real. Further, for zq > 1 and < zq < 1 Eqn. |4.4| 
returns values of a > 1 and a < 0, respectively. Consequently we need only consider real, 
negative values of the parameter z Q . 

We now turn to an examination of the saddle-point approximation of ( |3.9| ) for different 
regions of the phase diagram of the model at hand beginning with the weak coupling phase, 
H < 1. In this case the support of the eigenvalue distribution ( |2.9| ) is bounded away from 
9 = ±7r and hence the denominator in (|4.4[) is non-singular for all values of Zq. Consequently, 



in this regime a varies smoothly and monotonically with zq and the relation (|4.4j) can in 
principle be inverted to obtain z (a). With this information, the large N asymptotic form of 
the expectation value of the characters < Xa > can be determined by standard saddle-point 
methods. In Figure 3 we show a numerically calculated example of a as a function of zq for 
p = 0.5. For this same case we show a schematic diagram of the magnitude of the expectation 
value | < Xa > I as a function of a in Figure 4. In particular we see that the system has 
excitations in all irreducible representations. 

For p > 1 the situation is somewhat different. Now the support of the eigenvalue distri- 
bution (|2.9j ) is the full interval 6 G [— 7r, 7r], and the denominator of ( |4.4|) causes non-analytic 
behaviour to appear. As one increases p through unity the saddle-point relation for a shows 
this non-analytic behaviour as a discontinuity at z = — 1 (see Figure 5). The result is that 
an open interval of a values centered on a = 1/2 are mapped into this discontinuity when 
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Figure 3: Saddle-point relation for \i = 0.5 
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Figure 4: Schematic diagram of | < % Q > | vs. a for /x = 0.5 
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Figure 5: Saddle-point relation for \x = 1.2 
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Figure 6: Schematic diagram of | < Xa > \ vs. a for fi = 1.2 

the saddle-point relation ( |4.4|) is inverted. Since this discontinuity occurs in the saddle-point 
relation, it is not surprising to find that the curvature associated with the Gaussian integra- 
tion of the saddle-point approximation is divergent, effectively forcing the integral to vanish. 
In terms of the expectation values of different representations in the background of the non- 
Abelian Coulomb gas, we see that an open interval of fundamental representations centered 
about a = 1/2 is missing from the spectrum in the large N limit. In Figure 6 we show an 
example of the behaviour of the expectation value | < \ a > | with a for \i — 1.2. 

The main outcome of this analysis is that the expectation value of the central fundamental 
character < X1/2 > is vanishing if and only if fi > 1. Consequently it may be considered an 
order parameter distinguishing between the strong and weak coupling phases of the model. 
Physically the situation is clear. In the weak coupling phase the system can effectively screen 
the interactions of any pair of charges regardless of their representation since the system 
contains excitations in all representations of the gauge group. We conclude that the system 
looks much like a quark-gluon plasma where charges are effectively deconfined. At the phase 
transition line non-Abelian flux in the a = 1/2 fundamental representation becomes too 
energetically costly to produce and the system can no longer screen the interaction between a 
pair of a = 1/2 fundamental charges. In this strong coupling phase the interacting pair sees a 
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linear confining potential (though somewhat reduced as compared to the empty background). 
As one further increases \i the gap in the spectrum of fundamental representations becomes 
larger and in the extreme limit \x = oo the system contains only excitations in the trivial 
representation. This is precisely the confining phase of pure gluo-dynamics. 



5 Discussion 

As we have shown, the generalization of the concept of the Polyakov loop operator to probe 
the group theoretic excitations of a system of non-Abelian electric charges provides a conve- 
nient and unified way to quantify the physics of phase transitions. While the details of our 
presentation have centered on a two dimensional model with an infinite number of colours, 
the general concepts developed here should be applicable to interacting gauge systems in 
arbitrary dimensions for both infinite and finite rank (N) gauge groups. Unfortunately we 
are unable to test these ideas in a solvable two dimensional model as for finite N there is no 
phase transition and the system is always in the deconfined phase where all representations 
are present. 

One immediate problem with using the fundamental representations to characterize the 
phase diagram arises when considering finite, odd rank groups. For example in the physically 
relevant case of SU(3), there are only two fundamental representations k — 1, 2 and the 
order parameter < X1/2 > would naively denote the k = 3/2 fundamental representation. 
In terms of group theory this fractional representation is nonsense and strongly suggests 
that direct application of the large N results is not prudent. Even worse is the fact that a 
finite number of representations are not enough degrees of freedom to describe the change in 
character of the group elements from g ~ 1 in the weak coupling phase to g away from the 
identity in the strong coupling phase. Alternatively, we are free to use any independent set of 
irreducible representations to characterize the system. In particular the completely symmetric 
representations provide an equivalent algebraic basis to the fundamentals we have considered 
here. The strength of this approach is that there is no restriction on the number of symmetric 
representations for the unitary groups contrasting the N — 1 fundamental representations for 
SU(N). Unfortunately, repeating the calculations of Section 4 with symmetric representations 
one discovers that there is no strong signal for the phase transition and in particular one cannot 
define an order parameter. 

Despite the apparent difficulties with applying the current results to finite rank gauge 
theories, we feel that the general concept of characterizing the phases of a gauge theory 
coupled to matter by the spectrum of higher irreducible representations is useful. In particular 
it would be interesting to investigate these ideas in the setting of lattice calculations of gauge 
theories that are known to possess phase transitions. In fact, for the case of SU{2) gauge 
theory, the J = 1/2, 1 and 3/2 Polyakov loops have been calculated on the lattice |H| (see 



also 123). It would certainly be very instructive to have more complete information about the 
higher representation Polyakov loops in this simple model both with and without fundamental 
representation matter. 
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